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to introduce additional quadratic and linear time-dependnt terms to the given potential. 



S . In this way one can construct for general systems either coherent packets moving with 

^ . constant dispersion, or dynamically squeezed packets whose time-dependent dispersion 



> 

r> ' sponds to a suitable generalization of the displacement and scaling operators that generate 



remains bounded for all times. In the standard operatorial framework our scheme corre- 
sponds to a suitable generalization of the displacement and s( 
the coherent and squeezed states of the harmonic oscillator. 
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1. Introduction 

The present work addresses the problem of developing a comprehensive theoretical ap- 
proach to quantum control, in the sense of providing a unified dynamical principle to obtain 
the most general bounded quasi-classical wave-packet evolutions possible for arbitrary non 
harmonic and externally driven quantum systems. 

The developement of such a scheme is most easily available and conceptually clear 
by working in the framework of Nelson stochastic quantization. The latter is currently 
recognized as an independent and self-consistent formulation of nonrelativistic quantum 
mechanics in the language of stochastic processes [|l]-[H-[il. 

We remind that the experimental goal of quantum control is to use radiation to drive 
matter to a desired target or outcome [^, and several theoretical schemes modelling con- 
trolled wave-packet dynamics have been suggested. The unifying theme in all of the cur- 
rent theoretical and experimental work is the optimal use of the coherence of laser light to 
manipulate the quantum mechanical phase relationship among the eigenstates of matter. 

From this point of view, the coherent and squeezed states of the harmonic oscillator can 
be considered as special but paradigmatic examples of controlled wave-packet dynamics 
^ . Progress in the femtosecond pulse technology allows now for the realization of quantum 
control in the laboratory; for example, frequency-chirped femtosecond laser pulses have 
been synthesized to control the evolution of vibrational wave packets of the iodine molecule 

i- 

The potential interest of Nelson stochastic mechanics for the theory of controlled wave- 
packet dynamics stems then from the fact that the Nelson formulation of quantum me- 
chanics is nothing but a particular instance of classical stochastic control theory. Namely, 
it has been proven that quantum dynamics can be derived via a stochastic variational 
principle, by suitably extremizing the action functional of a classical mechanical system 
along diffusive trajectories replacing the classical deterministic ones. The stochastic vari- 
ational scheme has since been extended and exploited in a number of different contexts 
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For our purposes, what is appealing of the stochastic quantization scheme is the possi- 
bihty that it offers to obtain both new resuhs and new insights on old problems, by looking 
at quantum coherence in terms of general properties of classical diffusion processes. 

In particular, variational minimization of the stochastic osmotic uncertainty functional 
yields the complete structure of the quantum states of minimum uncertainty, i.e., the 
harmonic-oscillator coherent and squeezed states |l^. The same structure had been de- 



rived earlier |TT| by simply saturating the osmotic uncertainty inequality [|T^. 

In the present paper we then address, in the stochastic quantization approach, the 
problem of constructing a general framework for controlled wave-packet dynamics, pay- 
ing special attention to the construction of coherent and squeezed states for general non 
harmonic potentials. 

To that purpose, we introduce the conditions of classical motion for the wave-packet 
centre and the conditions of constant or bounded time-varying dispersion as constraints 
for the stochastic dynamics in an assigned, generic configurational external potential V{x). 
We then show that such constraints select a class of Nelson diffusions with classical current 



velocity and wave-like propagating osmotic velocity [0-[|14 

To each Nelson diffusion belonging to such class, there is associated a quantum state 
that is solution of the Schrodinger equation in a time-dependent potential V{x,t) having a 
simple relation with the original system V{x). In particular, if we ask for a purely coherent 
wave packet solution, i.e. whose centre (x) follows exactly the classical motion Xci{t) in 
V{x) with constant dispersion Ax, then the potential V{x,t) is completely determined by 
the functional form of V{x) and by the knowledge of Xci{t). 

Clearly, such "controlling" potential V{x,t) may be, in principle, experimentally fash- 
ioned in the laboratory, once a particular potential V{x) has been chosen. 

Connection with the standard operatorial language is then provided by showing that 
this new class of controlled coherent states is generated letting the standard unitary 



Glauber displacement operator act on any stationary state (for instance the ground state) 
associated to V{x). 

We then extend our scheme to construct coherent wave packets with bounded time- 
varying dispersion Ax. By resorting to the stochastic framework we are able to determine 
the new controlling potential V{x, t) (that in general does not coincide with the controlling 
potential V{x^t) obtained in the case of constant spreading) connected to the classical 
potential V(x), and to derive the evolution equation for Aa;, which turns out to be the 
classical envelope equation, well known in the theory of charged-particle beams dynamics. 

Knowledge of the solution Axd{t) of the classical envelope equation and of the classical 
trajectory Xci{t) associated to the potential V{x), determines unambigously the form of 
the controlling potential V{x,t). 

A suitable unitary operator acting on any stationary state (for instance the ground 
state) associated to V{x) is then introduced in the coordinate representation in order to con- 
struct these controlled states in the standard operatorial language: they are displacement- 
operator generalized coherent states with time-dependent dispersion. 

In fact, this unitary operator acts as the product of two distinct mappings on the 
ground state of the given potential V{x): the ordinary Glauber displacement operator and 
a dynamical scaling operator, namely a dynamical squeeze operator. Squeezing is then 
naturally embedded in this scheme, and the evolution equation for Ax yields also the 
dynamical equation controlling the time-evolution of squeezing. 

Overall, the above construction provides a natural way to introduce a class of physical 
coherent and squeezed states associated to general non harmonic potentials, in the sense 
of controlled wave-packet dynamics, by imposing suitable requirements on the physical 
properties of the desired solutions. In this way we provide a specific physical solution to 
the problem, posed by Schrodinger over seventy years ago, on how to generalize the notion 
of coherent state beyond harmonic systems. 

The paper is organized as follows. In section 2 we give a brief review of the basic 



ingredients of Nelson stochastic quantization that will be needed in the following. In 
section 3 we describe the structure of the harmonic-oscillator coherent and squeezed states 
in the stochastic picture. In section 4 we show how to construct general controlled coherent 
states in the stochastic framework by imposing the constrainsts of classical motion for the 
wave-packet centre and of constant wave-packet dispersion, and we derive the explicit 
expression for the controlling potential V{x,t). In section 5 we extend the discussion 
imposing the constraints of classical motion and of bounded time-varying dispersion and 
we derive the explicit expression for the displacement and squeeze operators associated to 
the controlled coherent states with bounded time-varying dispersion. In section 6 we draw 
our conclusions. 

2. Stochastic mechanics 

We shall quickly review the basic ingredients of the stochastic formulation of quantum 
mechanics that will be needed in the following. 

This quantization procedure rests on two basic prescriptions; the first one, kinematical, 
promotes the configuration of a classical system to a conservative diffusion process with 
diffusion coefficient equal to h/2m. 

If we denote by q{t) the configurational variable for a point particle with mass m, this 
prescription reads 



dq{t)=V(^+){q{t),t)dt+J—dw{t), dt>0. (1) 

In the above stochastic differential equation t;(+) is a (forward) drift ffeld that is determined 
by assigning the dynamics, and w is the standard Wiener process. 

An intuitive manner to look at Eq. (1) is to consider it as the appropriate quantum 
form of the classical kinematical prescription: the Wiener process models the quantum 
ffuctuations that are superimposed on the classical deterministic dynamics. 



If we consider the stochastic backward increment of the process dq{t) = q{t) — q{t — 

dt), under very general mathematical conditions the diffusion q[t) admits the backward 

representation 

dq(t) =V(_)(q(t),t)dt+\ dw*(t), dt > , (2) 

V 2m 

where w* is a time-reversed Wiener process. The probability density p{x, t) of the process, 
defined on the points x of the configuration space, is induced by the conditioned expecta- 
tions on the increments of the process with respect to the Wiener measure. The backward 
and forward drifts can then be expressed as stochastic fields on the configuration space 
through the conditioned expectations 



„,_,(x,t). JimJ«M-^|,(t)=x). (3) 

They represent respectively the mean forward (backward) velocity fields. 

In the hydrodynamic picture of the process, the drifts are replaced by the osmotic 
velocity u{x,t), 

u(.,t) ^ "'^''"■" ;"■-'<"■" = Av|,np(x.<)|^ (4) 

and by the current velocity v{x,t), 

v{x,t) ^ -i^^(-^^)+-i-)(-^^) _ (5) 

Finally, Fokker-Planck equation for the probability density p{x, t) takes the form of the 
continuity equation 

dtp{x,t) = -V[p{x,t)v{x,t)]. (6) 
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The process is completely determined by the the couple {v{x,t),u{x,t)), or, equivalently, 
by the couple {v{x,t), p{x,t)). 

The dynamical prescription is introduced by defining the mean regularized classical 
action A. In the hydrodynamic Eulerian picture it is a functional of the couple (p, f): 



A[p,v] 



ta 



m 



^{v\x,t)-u\x,t))-V{x,t) 



p{x,t)d^xdt , 



(7) 



where V{x,t) denotes the (possibly time-dependent) configurational external potential. 

The equations of motion are then obtained by extremizing A against smooth variations 
6p, 6v vanishing at the boundaries of integration, with the continuity equation taken as a 
constraint. 

After standard calculations one obtains the "quantum Newton law" 



dtv{x, t) + {v{x, t) ■ V)v{x, t) - -—V ; ' = -VV{x, t) 



Aw? 



\jp{x, t) 



with the current velocity fixed to be a gradient field at all points x where p{x,t) > 0: 



v{x,t) 



VS{x,t) 



m 



(9) 



where S{x, t) is a scalar field. 

Defining the wave function \l'(x, t) associated to a generic single-particle quantum state 
in the hydrodynamic form 



\E'(x, t) = J p{x,t)exp 



h 



S{x,t) 



(10) 



it immediately follows that the Schrodinger equation with potential V{x, t) for the complex 
wave function \l'(a;, t) is equivalent to the quantum Newton law together with the continuity 
equation, i.e. to two real nonlinearly coupled equations for the probability density p(a;, t) 



and for the phase S{x,t) (or, alternatively, for the osmotic and current velocities u{x,t) 
and v{x,t)). 

Then, to each quantum state there corresponds in stochastic mechanics a diffusion 
process q{t) with 

p{x,t) = \^{x,t)\\ (11) 



and 

^(x,t) 






In 



(12) 



\l'*(a;,t) 

The space integral of Eq. (8) yields the Hamilton-Jacobi-Madelung equation. It is 
useful for what follows to write this equation in the form 

dtS{x, t) + — t;^(x, t) u^{x, t) - -Vu{x, t) = -V{x, t) . (13) 

The correspondence between expectations and correlations defined in the stochastic 
and in the canonic pictures are 

{x) = E{q{t)), (p) = mE{v{x,t)) , 



Ax = Aq, (Apf = m^ [{Auf + [Avf] . (14) 

In the above relations x and p are the position and momentum operators in the 
Schrodinger picture, (■) are the expectations of the operators in the given state \E', E{-) is 
the expectation of the stochastic variables in the Nelson state {p,v} corresponding to the 
state \1', and A(-) are the root mean square deviations. 

In the theory of diffusion processes, the functional (Ag)^(A-u)^ is known as the osmotic 
uncertainty product; it shares the remarkable property that it is always greater or equal 
than the square of the diffusion coefficient. The following chain inequality then immediately 



follows: 



{Ax)\Apf > m\Aq)\Auf > —. (15) 



The osmotic uncertainty relation and its equivalence with the momentum-position un- 
certainty were proven in Ref . |T^] . 

3. Coherent and squeezed states of the harmonic oscillator 

Saturation of the osmotic uncertainty relation Eq. (15) yields the Glauber coherent 
states in the stochastic picture [|1T]: they are Nelson diffusions of minimum osmotic un- 



certainty, with constant dispersion Aq. They are characterized by a deterministic classical 
current velocity: 

v{x,t) = {v{x,t)) = Vd{t) = Xcl{t) , (16) 

and by an osmotic velocity linear in the process: 

In Eq. (16) and from now on we denote the stochastic expectations E{-) with the same 
symbol (■) used for quantum ones (wherever no confusion arises). We have also introduced 
the adimensional variable 

e = ^^^ , (18) 

obtained by first shifting the coordinate x by the classical trajectory Xci{t) = {q(t)), and 
then by scaling it through the wave-packet dispersion Aq. In fact, the adimensional 
configurational variable C, will play a fundamental role in all of the following. 

The probability density p{x,t) associated to the state (16)-(17) is readily obtained 
by comparing the expression for u{x,t) given in (17) and the relation (4) linking the 
probability density with the osmotic velocity. The phase S{x,t) is obtained comparing 
the expression for v{x,t) given in (16) and the relation (9) connecting the current velocity 

9 



with the phase. In this way one reconstructs completely the wave function \l'(x, t). The 
potential V{x,t), with classical trajectories Xdit), that solves Schrodinger equation for 
the state (16)-(17) is finally identified by replacing the expressions for S{x,t), v{x,t), 
and u{x,t) in the Hamilton-Jacobi-Madelung equation (13). This is a typical inverse- 
problems strategy: given a certain state with some desired features, in this case that of 
being of minimum osmotic uncertainty with constant dispersion, one looks for the potential 
that solves Schrodinger equation and allows for such state. 

We notice that the states of minimum osmotic uncertainty with constant Ag have 
a purely classical current velocity, so that Av = 0. Therefore the minimum osmotic 
uncertainty is exactly equivalent to the minimum Heisenberg uncertainty. Then, following 
the procedure outlined above, it is straightforward to recover the complete structure of 
the Heisenberg minimum uncertainty states, which are the Schrodinger-Glauber coherent 
states of the harmonic oscillator. We have: 



p(x, t) = , = exp 

/27r(Ag)2 



(X - Xd{t)f 

2(Ag)2 



Tfi 1 

S{x,t) = mvci{t)x - —Vci{t)xci{t) - -hut 



T / N 1 i (x — Xci{t)Y i 

m{x,t) = rexp<^-^ , /\V +T 

[27r(Ag)2]i "^ I 4(Ag)2 h 



f.-. ( Xd{t)\ h 



For those readers unfamiliar with the stochastic picture, we remind that in the above 
expressions (19) Xci{t) = (x) and mvci{t) = (p). 

If we now impose saturation of the osmotic uncertainty, but allowing for a time- 
dependent dispersion Ag, we obtain the harmonic-oscillator squeezed states. They are 
Nelson diffusions with time-varying Ag, with osmotic velocity u{x,t) still of the form 
(17), and with current velocity allowing for a term dependent on the time-evolution of the 
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dispersion |n]]: 

v{x,t) = v,i{t) + ^-Aq . (20) 

at 

The last term in Eq. (20) is responsible for the quantum anticommutator term ap- 
pearing in the phase of the squeezed wave packets, which are quantum states of minimum 
Schrodinger uncertainty. This can be easily seen as follows: define the centered position 
operator x^ = ^ — {x), the centered momentum operator pc = P — {P) , and the centered 
Nelson process qdt) = q{t) — {q{t)). All this quantities have zero expectation value. Next, 
consider the quantity qcit); obviously, {q^(t)) = (Ag)^. By straightforward calculations 
one gets 

^{q'M = ^(A?)' = 2[(g(t)t;(a:,t)) - {q{t)){v{x,t))] = ii^^iM. , (21) 
at at m 

where {, } denotes the quantum anticommutator, which expresses, when taken between x^ 
and pci the quantum position-momentum correlation. In turn, the latter is the Schrodinger 
part of the uncertainty, and for the squeezed states it is directly connected with the cur- 
rent uncertainty product in the Nelson stochastic picture. In fact, from Eq. (20) it is 
straightforward to see that 

{Avf = 2[j^Aq^ . (22) 

The classical evolution equation for the dispersion Ag can be easily obtained by inserting 
the expressions (17) and (20) for the stochastic hydrodynamic velocities into the quantum 
Newton law (8) and the Hamilton-Jacobi-Madelung equation (13); after straightforward 
manipulations, we obtain 



d'Aq 



+ muj'Aq = , (23) 



df^ ^ 4m2(Ag)3 ' 

which is the classical equation for the beam envelope, well known from the theory of 
classical optics and of particle accelerator dynamics. 
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The expressions (19) for the Gaussian probabihty density and for the harmonic oscillator 
potential remain unaltered in the case of minimum uncertainty states with time-dependent 
dispersion. However, from Eqns. (20) and (21) it follows that the phase of the wave function 
picks up an extra term proportional to the quantum anticommutator: 

'if(x,t) = r exp 

[27r(Ag)2]3 

f {q{t)v{x,t)) - Xci{t)vci{t) \ ^ ^c^K-j ... . . 

m\ — ) — - —Vd{t)xd{t) - -ut ). (24) 

The above stochastic picture for the harmonic-oscillator coherent and squeezed states 
can be entirely derived in a stochastic variational approach by extremizing the osmotic un- 
certainty product against smooth variations of the density p{x, t) and of the current velocity 
f (x, t) ||10|. The possibility of extending such variational approach to study local minimum 



{X-Xd{t) 


'-^, 


'm 

—Vcl{t)X + 


I 4(Ag)2 


X - Xd{t) 

[ Aq \ 


2 


m 

Vrl 


it)xd{t) - \ujt 



uncertainty behaviours in non harmonic systems is a current subject of investigation |T0 

This concludes our discussion of the coherent and squeezed states of the harmonic 
oscillator in the framework of stochastic mechanics. We next move to study the problem 
of how to generalize the concept of coherent and squeezed states to non harmonic systems, 
from the same perspective of controlled wave-packet dynamics that we have adopted in 
this section, with the emphasis focused on the determination of the potential that solves 
Schrodinger dynamics once a certain state with some desired features has been selected. 



4. Controlled coherent quantum wave packets: constant dispersion 

A well known property of the coherent states (16)-(17) and of the squeezed states 
(17)-(20), is that they follow the classical motion in the coherent Glauber sense: 

j^{m{v{x,t))) = -VV(x,t)U=(g(i)) . (25) 
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Again, from a dynamical point of view a coherent state is a wave packet whose cen- 
ter follows classical motion not only in the mean, but even along the mean (classical) 
trajectories, and whose dispersion is either constant, or controlled in its time-evolution 
(squeezing). 

In quantum mechanics the dynamics of mean values obeys Ehrenfest theorem: as a 
consequence, the coherent evolution Eq. (25) is satisfied exactly if 

{VV{x,t)) = VV{x,t)U=^,) . (26) 

In the case of quadratic potentials the above constraint is authomatically satisfied for 
any quantum state. For other generic potentials V{x,t) Eq. (26) in general cannot be 
satisfied. Our strategy is now to search, in the framework of Nelson stochastic mechanics, 
for states obeying the coherence constraints of classical motion (25)-(26) with constant 
dispersion Ag. 

To this purpose, we first observe that Eq. (21) can be recast in the form 

Ag^ = {q.{t)-v{x,t)), (27) 

where the dot denotes the scalar product between the centered process qdt) and the current 
velocity v{x,t). 

Therefore, in stochastic mechanics, we have that necessary and sufficient condition for 
a constant dispersion Ag is that the expectation value of the scalar product between the 
centered configurational process and the current velocity vanishes. We see immediately that 
a sufficient condition for this to happen is that the current velocity be purely classical, that 
is, v{x,t) = Vciit). Other possible choices of current velocities orthogonal to the centered 
Nelson process could be in principle considered, and they might lead to the definition of 
new classes of states in quantum mechanics. For the moment being, we concentrate on the 
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simplest choice of a classical current velocity. 

In stochastic mechanics, due to the existence of the continuity equation, a particular 
choice of the current velocity selects an entire class of osmotic velocities, and thus of 
quantum states. In particular, the choice v = Vci{t) that guarantees a constant dispersion, 
does not restrict the osmotic velocity to be of the minimum uncertainty form (17). Rather, 
upon insertion of Vd (t) in the continuity equation, one can show that the latter is satisfied 
by any osmotic velocity (probability density) of the following wave-like propagating form 



where G{C,) can be any arbitray adimensional function of the adimensional shifted and 
scaled configurational coordinate ^, provided it yields a normalizable probability density 
p(^). Keeping in mind relation (4) connecting the density to the osmotic velocity, and re- 
minding that a probability density (being non negative) can be expressed as the exponential 
of a real function, we have: general form 

P(0 = ^exp[2i?(0]. (29) 

In the above. A/" is a positive, adimensional normalization constant, and d denotes the 
spatial dimension of the system under consideration. Being V^. = V^/Ag, it follows that 
the adimensional normalizable function R{C,) is related to G{C,) by G{C,) = 2V^-R(^). 

The associated wave function is thus of the coherent-state form 



ArV2 



R{0+^(jvAt)-x+^' 



(30) 



where we remind that mvci{t) = (p) and 5*0 (t) is an arbitrary time-dependent constant. 
What are the properties of this class of states selected by the couples (f , u) of the form (16)- 
(28)? By construction, they are nonstationary states of constant dispersion. Of course, in 
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general they are not Heisenberg minimum uncertainty states; the latter are recovered only 
with the choice of the linear form G{C,) = — ^. 

However, it is easily verified that they share with the harmonic-oscillator minimum 
uncertainty states the remarkable property that they obey the constraint (25) for classical 
motion, apart, at most, a constant. To see this in the simplest and most explicit way, let 
us suppose that Vdit) is chosen such that its integrals Xdit) are the classical trajectories 
of a generic one-dimensional configurational potential V{x). What follows can then be 
generalized with minor technical complications to higher spatial dimensions. 

Suppose then that we are given the stationary states of the quantum mechanical prob- 
lem associated to V{x). For instance, we consider the ground state \l'o(a;,t). By a simple 
scaling argument, it can be cast in the general form 



vAgo 



X \ I ^ 



(31) 



where A/o is the adimensional normalization constant. Ago the ground-state dispersion, Eq 
the ground-state energy, and F a given adimensional function. 

Comparing expressions (30) and (31) we see that the coherent wave function '^d^^t) 
is obtained from the ground state wave function '^^{x^t) first by identifying the arbitrary 
function R with F, and then by applying to \l/o the following unitary Glauber-like dis- 



placement operator [M: 



I)[xd{t),Vd{t)] = exp 



USo{t)-Eot) 
n 



exp ( -mvci{t)x] exp {--Xd{t)p ) . (32) 



This operator, when applied to any wave function \l'(x,t) displaces its space argument 
X into X — Xdit) and introduces in the phase the coherent term mvd(t)x + S'o(t) — E^t. 

We now see that the wave packet (30) is really a coherent state in the sense that: the 
associated probability density has the same functional form of the ground-state density. 
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so that it shares the same statistics; in particular, the two normahzation constants co- 
incide and the constant spreading Ag is just the constant ground-state dispersion Ago- 
Furthermore, the uncertainty product remains constant too and equal to the ground-state 
uncertainty product, just like the standard harmonic-oscillator coherent states do; in fact 
the latter are just a particular case of the present construction. Last but not least, the 
wave-packet center follows the classical motion Xci{t) in the given configurational potential 
V{x) according to the Glauber law (25). 

The price to be paid for this construction is that these states do not satisfy the time- 
dependent Schrodinger equation in the originally assigned potential V{x), unless of course 
the latter is chosen to be a polynomial of degree not greater than two. 

However, the controlled states "^dxjt) of the coherent form (30) are still solutions of 
the time-dependent Schrodinger equation in a modified potential V{x,t) which has a very 
remarkable relation to the original potential V{x). Namely, taking the wave-like density 
p{^) Eq. (29) and the coherent phase Sc{x, t) = mVci(t)x+So(t) associated to the state (30), 
and inserting them in the Hamilton-Jacobi-Madelung equation (13) one finds that "^d^, t) 
is a solution of the time-dependent Schrodinger equation in the following time-dependent 
potential: 

V{x^t) = V[x — Xci(t)] + mxci(t)x . (33) 

The above expression gives the generic form of a controlling potential that allows for 
the desired wave-packet solution with coherent and localized dynamics in some previously 
assigned configurational potential V{x). 

Formally, the controlling potential V{x,t) is obtained from the original external poten- 
tial V{x) shifting its argument by the classical trajectories associated to V{x) itself, and 
then by adding a correcting term linear in x, which is multiplied by the inertia associ- 
ated to the force field —VxVlx). The time-dependence of the controlling potential is thus 
parametric via the solutions of the classical equations of motion in the external field V{x). 
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The question naturally arises about the physical interpretation of this scheme of quan- 
tum control. In particular, one may ask what is the meaning of the correcting time- 
dependent linear term and more generally whether a controlling potential of the form (33) 
might be fashionable in the laboratory. 

In principle, the potential (33) is a well defined object, and the practical realization 
of it heavily depends on the actual choice of the original potential V{x): one should then 
take up a careful case-by-case analysis. However, we can give some general comments on 
the structure of V{x,t) for arbitrary V{x). 

We observe that the linear correcting term superimposed on the original potential, has 
the simple interpretation of an electric field with amplitude varying in time according to 
the classical force law in V{x). The problem is then what actual time-dependences can be 
experimentally realized. As for the shift in the argument of V{x), it is in principle feasible, 
provided again that the assigned time-dependence of the coefficients multiplying powers 
of X can be actually fashioned by some wave-form generating set up; yet, we observe that 
there might be many situations where the shift could be considered a small perturbation, 
allowing for Taylor expansion of V or other approximate treatements, for instance in the 
case of a dynamics taking place in spatial and/or time regions where Xci{t) is a slowly 
varying quantity. 

In conclusion, the states that we have constructed in the stochastic framework are 
controlled wave packets with optimized quasi-classical behavior, that generalize the con- 
cept of the harmonic-oscillator coherent states in a precise physical sense. They follow 
a classical motion with constant dispersion in a given configurational potential V{x) and 
can be obtained via a unitary transformation of a generic energy eigenstate of V{x). The 
programming potential V{x,t) that must be introduced for the desired states to satisfy 
quantum dynamics has a simple and intriguing structure, strongly related to the original 
potential V{x). In particular, its form allows in principle for approximate treatements 
according to what degree of localization and coherence one desires to accomplish for a 
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generic nonharmonic quantum system. 

5. Controlled coherent quantum wave packets: time— dependent dispersion 

We now proceed to consider the case of time-dependent dispersion, that is we consider 
the more general form Eq. (20) for the current velocity of minimum uncertainty, and we 
follow again the strategy adopted in the previous section for the case of the classical current 
velocity with constant dispersion. 

We expect that the controlled coherent states that we will select by taking the choice 
(20) for the current velocity should be states following two coupled dynamical equations, 
Eq. (25) for the wave packet center, and an evolution equation for the dispersion Ag, 
analogous of the envelope equation (23) derived in the harmonic case. 

We proceed as follows. We first set the notation by relabelling the dispersion as a 
function of time. We put: 

a{t) = Aq . 

Next, inserting Eq. (20) into the continuity equation (6) we are left with: 

aM..t) = "«■') W) - -it) ^ (34) 

a{t) 

whose general solution is function only of C,, as can be immediately seen, e. g. by moving 
to Fourier space. Again, one has selected a class of probability densities of the form (29) 
as well as a class of osmotic velocities of the wave propagating form (28). 

Inserting now the current velocity (20) in the equation of motion (8), by Eq. (28) we 
obtain 

- mia{t) + -V.u\i) + -Vlu{0 = V.\/(x, t) - (V.\/(x, t)) , (35) 

where we exploited Ehrenfest theorem Vci{t) = —{'VxV{x,t))/m. 
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Letting x = {q{t)) = Xci{t) (i.e. ^ = 0), we are left with 

V,\/(x,t) U=,, -(V,\/(x,t)) = -V^u\i) \^^,+-Vlu{i) 1^=0 ■ (36) 

It is straightforward to show that this relation holds also in the case of constant dis- 
persion v{x,t) = Vci{t). The right hand side is obviously either constant or zero except 
for singular potentials: in these cases u{C,) diverges in ^ = 0. However, the scheme can 
be implemented also for singular potentials by taking ^ = Xci{t) rather than x = Xci{t). 
Explicit examples and applications to both singular and non singular potentials will be 
discussed in detail elsewhere. 

Given the state (20)-(28) we now want to write explicitely the associated phase S{x, t) 
and the evolution equation for the dispersion 0"(t). This is achieved by exploiting Hamilton- 
Jacobi-Madelung equation, Eq. (9): reminding that v{x^t) is the gradient field of S{x,t), 
Eq. (20) implies 

S{x,t) = mv,i{t)x + !lfc-^lM!^(t) + So{t) . (37) 

Inserting Eqs. (37), (20), and (28) in Eq. (13), and taking its expectation, we obtain: 

: -{V{x,t))+x,i{t){V^V{x,t)) - Soit) . 



-a{t)ait) + - 



^2 ^/ 2, 



{v{u)r - j{u\o) 



We can eliminate in the above expression the time-derivative of the classical action and 
obtain, after trivial manipulations, an evolution equation for o"(t): 

,(t)-W««=-«Z£lW)). (38) 
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By Eq. (28) for u{^) it is immediately seen that 



nm 






(39) 



where K"^ = /_^ G^iOpiO^^i Eq- (38) is then the desired equation for the time-evolution 
of the dispersion, generalizing the classical envelope equation (23). Moreover, it is coupled, 
through the gradient of the potential V{x, t), with the equation of motion (25) for the wave 
packet center {q(t)) = Xdit). 

The general form of the wave function for this class of states is readily obtained putting 
together Eq. (29) for p(0 and Eq. (37) for S{x,t): 






a 



(^)]d/2 



RiO + I {rnvd{t)x + ^^^I^^^^y^^W + S,{t) 



(40) 



We can rewrite this expression in more standard quantum mechanical terms by remind- 
ing Eq. (21), so that 



"^six.t) = 7777^7777^ exp 



o- 



(t)]^/2 



^(0 



h 



{p)x 



{{QcPc}) 

{2a{t)f 



x-{q)f + Soit) 



(41) 



The above wave functions "^sixjt) are realizations of nonstationary states with classical 
motion and controlled time-dependent spreading. They generalize the harmonic-oscillator 
squeezed states (24) in the same sense as the states (30) \E'c(a;, t) generalize the harmonic- 
oscillator coherent states (19). 

We shall now move to the study of the controlling potential that needs to be introduced 
for the controlled squeezed states \E's(a;, t) to satisfy Schrodinger equation. But first we wish 
to complete the analogy with the controlled coherent states introduced in the previous 
section. Namely, we will show that the states (40) can be obtained introducing a proper 
scaling "squeeze" operator. 

We proceed as follows. We first recall that the harmonic-oscillator squeezed states are 
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generated by the successive application of a scaling squeeze operator and of Glauber dis- 
placement operator to the harmonic-oscillator ground state. We then define the following 
dynamical scaling operator 5'[cr(t)]: 



5'[cr(t)] = exp 



Hfa« 



exp 



'im 



exp 



/W^Wrr. .T .21 



2n{aoy 



■[{q,p},Q 



(42) 



where uo denotes the (time-independent) dispersion associated to the ground state "^oi^, t) 
of a certain assigned configurational potential V{x). Given a{t) solution of Eq. (38), the 
two functions f{t) and g{t) read 



fit) 



2 I ao , 



771 n 



(43) 



We see from these relations that the function f{t) plays the role of a dynamical squeezing 
parameter. We now let 5'[cr(t)] act on the ground state wave function \l/o(x, t), cast in 
the form (31), associated to a given configurational potential V{x). We so define the 
dynamically scaled wave function 



^sc{x,t) = S[a(t)]-^o(x,t). 



(44) 



By straightforward algebra,{g,p} = ih{l + 2xd/dx) and [{^,p},^^] = — 4i/lg^, and one 
obtains 

X{x,t), (45) 



^!sc{x,t) = exp 
where xi^^t) is given by 



fit) (l + 2x^ 



X{x,t) = exp 



#(1 - 2/w)-^' 



^o(a;,t) 



(46) 



We now exploit the extension to the real axis of the following relation, holding for any 



analytic function W{z), that was introduced in Ref. [^ for the study of g-oscillators 
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coherent states: 

Q-± [W{x)] = W{Qx) , (47) 

with Q a real c-number and W analytic on the real axis. Letting Q = exp[2/(t)] and 
W = X one finally is left with 



'^sc{x,t) = exp 



/(t)+.^(l-2/(t))e^/(V 



^o{e^f^'^x,t) . (48) 



We then obtain the state \E's(a;,t), Eq. (40), by applying D[xci{t),Vci{t)], Eq. (30), to 
^,,(x,t): 

^,(x, t) = D[x,i{t),v,i{t)] {S[a{t)] ■ ^o{x, t)) = D[x,i{t),v,i{t)] ■ ^,,(x, t) . (49) 

By recalling Eqs. (43) it is then straightforward to show that ^^(a;, t), Eq. (49), 
coincides with the wave function (40). We thus proved that the controlled squeezed states 
(40)-(41) can also be introduced by a suitable modification of the displacement-operator 
approach to the harmonic-oscillator coherent and squeezed states. 

The states \E's(x, t) do not satisfy Schrodinger equation in the assigned potential V{x), 
although their wave packet center does follow the classical motion in the force field gener- 
ated by V{x). 

The situation is completely analogous to that described in the previous section. How- 
ever now, due to the presence of a time-varying dispersion, the controlling potential does 
not coincide with the expression (33). We thus label it V{x, t) to distinguish it from V{x, t) 
introduced in the previous section. After solving the Hamilton- Jacobi-Madelung equation, 
the new controlling potential reads: 

V{x,t) = V[x - x,i(t)] + m (i,i{t) - ^x^iit)] X + ^^^a;2 . (50) 
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We see that, compared to the controlled coherent case, a controlled squeezed state must 
be associated to a programming potential with an extra correcting quadratic term; further- 
more, the time-dependent coefficients of the correcting terms acquire a more complicated 
structure as they now depend not only on the classical trajectories Xci{t) but also on the 
solutions a(t) of the generalized classical envelope equation (38). 

The physical interpretation of the quadratic correcting term can be simply given in 
terms of a diamagnetic interaction superimposed by the controller in the laboratory to the 
previously existing external interaction V{x). About the experimental feasibility of such a 
controlling set up we can repeat in principle what we have observed in the previous section 
for the linear controlling potential. We will further comment on both the linear and the 
quadratic controlling potentials in our conclusions. 

6. Conclusions and outlook 

Let us summarize at this point our results. 

Working in the framework of Nelson stochastic quantization, we have introduced a class 
of controlled coherent states (constant dispersion) and a class of controlled squeezed states 
(bounded time-evolution of the dispersion). 

The wave-packets centers follow a classical evolution in a generic preassigned time- 
independent external potential V, while the states obey Schrodinger dynamics in a time- 
dependent potential that has a simple relation to V. For constant dispersion such time- 
dependent potential (that we have named controlling or programming potential) presents 
a linear correcting term to the original time-independent potential, see Eq. (33). For 
the bounded time-dependent dispersion, the controlling potential presents a linear plus a 
quadratic correcting term, see Eq. (50). 

For the controlled squeezed states, the evolution equation controlling the spreading of 
the wave packet is the classical envelope equation of classical optics, and it is naturally 
coupled with the classical evolution equation for the wave-packet center. The solutions 
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of the two classical equations enter as (time-dependent) parameters in the programming 
potentials. 

We have also showed that the controlled states can be obtained through a particu- 
lar extension of the displacement-operator coherent and squeezed states of the harmonic 
oscillator by defining a suitable dynamical scaling (squeeze) operator. 

The consequent dilatations or contractions of the wave-packet width are then shown 
to be controlled by a single adimensional squeezing parameter f{t). 

A comment is due at this point. It is well known that generalized coherent states 
can be obtained extending the three different existing approaches to the definition of the 
harmonic-oscillator coherent states: they are, respectively, the minimum-uncertainty, the 
annihilation-operator, and the displacement-operator method. 



The states obtained by extension of these three methods are in general different [T^ . 

The displacement-operator coherent states are those preserving most of the properties 
of the harmonic-oscillator coherent states: they are still overcomplete and still enjoy res- 
olution of unity. Moreover, we have shown in the present paper that by a proper choice 
of the Glauber parameter a entering the Glauber displacement operator, they follow a 
classical motion without dispersion. 

As to squeezed states, an extension of the minimum-uncertainty and annihilation- 
operator methods to arbitrary nonlinear systems was carried out by Nieto and collaborators 



|T6| . They also introduced an extension of the minimum-uncertainty method to obtain 



generalized coherent states 17 



However, an extension of the displacement-operator method to obtain generalized 
squeezed states runs into difficulties [|18| and is still missing, although some progress in 
that direction has been recently obtained by Nieto and Truax for systems allowing for 
Holstein-Primakoff or Bogoliubov transformations ||T9 |. 



We have here introduced a possible extension of the displacement-operator method 
by constructing squeezed states via a structure of controlling potentials. In this sense, 
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we have given a particular answer to the problem posed by Schrodinger exactly 70 years 
ago, whether it is possible to construct coherent-states solutions for general non quadratic 
potentials. We showed that it is possible, at the price of modifying the original potential 
in a very definite way: through the external action of a correcting interaction, the original 
system can follow exactly and indefinetely a coherent (or a squeezed) dynamics. 

In this paper we outlined the general features of the method: elsewhere [^ we will 
present the construction of coherent and squeezed states in the sense of Schrodinger for a 
wide sample of potentials of physical and conceptual interest. 

It might be worth noting that we have chosen for simplicity the ground state to gen- 
erate controlled coherent and squeezed states by the action of operators (32) and (42); it 
is however immediately seen that their application on any stationary state yields again 
controlled coherent and squeezed states of the form (30) and (40). 

We remark that Eq. (39) represents the "stochastic squeezing" condition satisfied by 
our states. Namely, it expresses the complementary time-dependence of the spreading a(t) 
and of the osmotic velocity uncertainty Am. 

It is easily seen that in the canonic picture Equations (14), (20), and (39) imply 
Ag2Ap2 = K'^ + L'^(t), with L{t) = mAqd{Aq)/dt. 

The reciprocal variation in time of Ag and Ap is then ruled by Ag itself, determined as 
the solution of Eq. (22) with the initial condition Ago. In this way squeezing is introduced 
as a self-consistent prescription on the dynamical evolution of the wave-packet spreading. 

In conclusion, our scheme of coherent and squeezed states via programming potentials 
provides an instance of the so-called theory of quantum control (or "Controlled Quantum 
Mechanics"), in the sense that, given a desired quantum solution (e.g. a squeezed state), 
one can provide a theoretical framework to describe what dynamical system (potential) 
must be introduced to produce such a state. 

This theoretical model is deeply rooted in the ideas and techniques of the theory of 
stochastic optimal control: a probabilistic approach to the description and the construc- 
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tion of quasi-deterministic structures, as generalized coherent and squeezed states, must 
involve an optimization procedure (minimization of the noise) and an external dynam- 
ical monitoring (programming potentials) of the non-deterministic system under study 
(holding the wave packet localized along the classical trajectory). 

We are currently studying the possibility to apply our scheme to two very actual and 
important problems in the field of quantum control: Rydberg wave packets ||22|-||^ and 



molecular pseudo-gaussian states |^ . In both instances, the experimental situation consists 
of a laboratory-fashioned interaction (femtosecond laser pulses) superimposed on existing 
natural interactions (the atomic and molecular potentials) that gives rise to an approxi- 
mately coherent and localized dynamics of the wave packets. 

Our scheme of quantum control seems potentially able to suggest from a general the- 
oretical framework the optimal interactions that should be fashioned to obtain the best 
control on the dynamics (the highest degree of coherence and localization). Namely, it 
would be interesting to verify the actual experimental feasibility of the programming po- 
tentials (33)-(50) and of the associated controlled coherent and squeezed states for some 
specific physical systems such as the ones mentioned above. We will report elsewhere [73 
about the work currently in progress on these subjects. 
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